Relationship between dynamical heterogeneities and stretched exponential relaxation 
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We identify the dynamical heterogeneities as an essential prerequisite for stretched exponential 
relaxation in dynamically frustrated systems. This heterogeneity takes the form of ordered domains 
of finite but diverging lifetime for particles in atomic or molecular systems, or spin states in magnetic 
materials. At the onset of the dynamical heterogeneity, the distribution of time intervals spent in 
such domains or traps becomes stretched exponential at long time. We rigorously show that once 
this is the case, the autocorrelation function of the renewal process formed by these time intervals 
is also stretched exponential at long time. 
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I. INTRODUCTION 

The stretched exponential relaxation (SER), I(t) ~ 
cxp[— (t/r) 13 ], describes the relaxation of pure glasses 
very accurately, often over very wide ranges [IJ. For ex- 
ample, SER describes equally well stress experiments on 
amorphous Se, centered on 10 3 s, as well as spin-polarized 
neutron scattering experiments, centered on 10 -9 s, both 
with the same value of /3, over a time range of more than 
10 12 . SER also describes accurately stress relaxation 
of the superconducting transition temperatures in the 
cuprates, probably associated with the two-gap domain 
structure revealed by high resolution scanning tunneling 
microscopy |2j; moreover, the measured value of (3 is the 
predicted value, 3/5, for YBCO, by far the most exten- 
sively studied cuprate, with the highest quality samples. 
The agreement is such that it has been recently proposed 
that the SER be used as an independent measure of sam- 
ple quality in these materials Q. Phillips' review pj gives 
many more examples and demonstrates the power of SER 
as a guide to understanding the dynamical properties of 
intrinsic glasses, especially network glasses. One of the 
most spectacular examples was the successful prediction 
P| of the long- and short-preparation-time values of (3 of 
orthoterphenyl (OTP), as measured by multidimensional 
NMR. OTP is the purest organic glass former available 
commercially; the predicted values were 0.43 and 0.60, 
and the observed values were 0.42 and 0.59. 

It appears that SER may arise as a result of spatial 
inhomogeneities that are quenched into macroscopically 
pure, non-phase separated glasses on a microscopic scale 
as traps or sinks that appear as part of the glass-forming 
process. The existence of these inhomogeneities is dif- 
ficult to detect experimentally, but they readily account 
for SER, and predict correctly the values of beta in many 
experiments 1]. Dynamical heterogeneities have been 
seen experimentally and numerically in glasses and su- 
percooled liquids 0,13, spin-glass models yllM, Lennard- 
Jones alloy mixtures with ellipsoidal probes |9j as well as 
in coupled-chaotic systems such as diode-resonator ar- 
rays 10]. In particular, the concept of dynamical het- 



erogeneities has emerged as critical for the description of 
the microscopic dynamics associated with the dramatic 
slowing down of the relaxation and diffusion processes 
as the temperature approaches T g , the glass-transition 
temperature [Tl 111 . 

These spatial heterogeneities occur over time scales 
that are long with respect to the basic unit time of the 
systems — one step in coupled maps and a lattice vi- 
bration in supercooled liquids — but shorter than the 
macroscopic relaxation, where materials are known to be 
homogeneous. As the temperature or the driving force 
is decreased, these regions become more and more stable 
and might overtake the whole sample at a transition tem- 
perature. Although the nature of spatial heterogeneities 
and their relation to the macroscopic dynamics remains 
incompletely understood, there is a considerable amount 
of work that relates directly to this issue. The mode- 
coupling theory describes the dynamics from the liquid 
phase into the supercooled regime and makes quantita- 
tive predictions that have been verified numerically in 
many systems [3 0|. It cannot yet describe accu- 
rately the very long time behavior in the supercooled 
regime, however |l4j. Numerically, dynamical hetero- 
geneities have been studied using a wide r ang e of criteria 
tied various aspects of this feature^, lla fl7| . Recently, 
however, some groups have focused on the renewal or 
mean first-passage time (MFPT) and the waiting- 
time distribution (WTD) which provide some link to the 
structure of the energy landscape 0, |2(J . 

Although both the dynamical heterogeneities and the 
stretched exponential behavior are characteristics of 
many frustrated systems (e.g. supercooled liquids |2l|). 
it is not clear what the relation between these two prop- 
erties is. In this paper, we address this question and show 
that dynamical heterogeneities can be directly responsi- 
ble for the observation of stretched exponential relaxation 
in coupled map lattices. We also show that these results 
are applicable to supercooled liquids and can provide a 
microscopic basis for the stretched exponential relaxation 
in these systems, in line with the observations of Denny et 
al. on the long time dynamics of meta-basin hopping |20| . 
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Recently, Hunt et al. found that the distribution of the 
renewal time of a certain stochastic process measured in a 
one-dimensional lattice of coupled diode resonators could 
be fitted to a stretched exponential function over 6 or- 
ders of magnitude [Toj |. Simulations on a related model 
of diffusively coupled via short range interactions non- 
linear maps show that the fit can be extended to more 
than 9 orders of magnitude, ruling out any other power 
law or simple combination of pure exponentials [lfl [23 • 
The quality of the fit in this system is sufficient to dis- 
tinguish the region of the parameter space where the dy- 
namics is described by power-law distribution coupled 
with an exponential cut-off from the really stretched ex- 
ponential distributions. (Ref. [^ contains background 
information important for understanding the results of 
the present work and will be referred to as [I] in the fol- 
lowing.) 

These simulations [I] also demonstrated the impor- 
tance of the system's size on its dynamics. While 
stretched exponentials are sufficiently robust and are ob- 
served even in some small samples, the quality of such 
fits is much poorer than in the big samples. This is in 
direct analogy with real experimental results and sup- 
ports the observation that the quality of the samples is 
crucial for obtaining unambiguous stretched exponential 
behavior and meaningful values of the exponent f3 0, Q . 

Using simulations on this and other models, as well 
as some analytical results, we show here that (1) the 
stochastic process studied in Refs. ^(j and [I] can be 
identified with a renewal process [2^, (2) a stretched ex- 
ponential distribution of the renewal time implies a sim- 
ilar shape for the decay of the two-point auto-correlation 
function, and (3) renewal processes are present in a su- 
percooled liquid and play a central role in its dynamics. 

Granted the universal features of SER, readers unfa- 
miliar with chaos theory and coupled maps may not see 
immediately the relevance of these mathematical tools to 
understanding the dynamical properties of supercooled 
liquids and glasses. Even the successful predictions of the 
trap (sink) model 0, |3| appear to be unexpected. For 
instance, the survival probability of a random walker in 
the presence of randomly distributed static traps crosses 
over from exponential to SER only at longer times where 
very few (~ 10~ 30 for 10% traps) of the walkers have 
survived |25|. In our examples of steady-state dynamics 
in coupled maps, the analogous crossover in trap-time 
distributions p(t) occurs near KT 1 to 1CT 2 x max{p(i)} 
depending on the value of a control parameter. Even 
the latter number is small compared to the crossovers 
observed, for example, for relaxation of the first peak in 
the structure factor in molecular dynamics simulations 
(MDS) of crystallization-avoiding soft sphere mixtures 
[26| . which occur near 3/4 of the maximum value. The 
reason for these differences is that the systems chosen 
for MDS already represent prepared states with strong 
glass-forming tendencies. The random walker with ran- 
dom traps model contains spatial inhomogeneities, which 
we will show represent an essential feature necessary for 



SER, and which are absent from some popular models of 
static glass structure such as mode-coupling theory [27| . 
That model does not contain the collapse of phase space 
that occurs as the liquid is supercooled to the glass tran- 
sition, and that is why SER appears only at very long 
times. However, it does contain the correct dimensional 
dependence of diffusive behaviour that leads to the rela- 
tion (3 = d/(d + 2), which is important for comparisons 
to experiment. That the one-dimensional coupled map 
lattice discussed here reaches a steady-state much more 
rapidly than the random walker is a clear indication that 
this model is physically relevant. The way in which this 
happens was illustrated in Fig. 8 of [I]: the surviving 
periodic domains are those that have avoided chaotic re- 
gions. 

We also emphasize that in this simple model, the be- 
havior of the stretching exponent j3 as a function of a 
control parameter r is similar to the dependence of (3 on 
the temperature in both 3D [28j and ID [8j spin-glass 
models. This similarity not only captures the qualitative 
tendency for (3 to decrease with decreasing r (or T), but 
also the quantitative agreement in the low-r (or T) limit 
(3 w 1/3. The same lowest value of_B was obtained in 
studies of random walks on fractals [23. l30l l3l| . 

II. RESULTS 
A. Properties of the coupled map model 

The nonlinear model of Ref. ^(j; used to repro- 
duce qualitatively the experimental system, is a one- 
dimensional chain of N diffusively coupled nonlinear de- 
terministic maps, f(x), with a coupling constant a and 
periodic boundary conditions. The time evolution of this 
system is discrete and is described by the following iter- 
ative equation, 

X n (t+1) = (l-a)f{x n (t)) + ^ {/(z B -l(*)) + f(Xn+l(t))} ■ 

(1) 

which, given initial conditions x n (0) for each site n = 
1, 2, . . . , N, generates a time series {x n (t)}. The inter- 
action is totalistic and involves only the nearest neigh- 
bors. Following [I], in this study we use the logistic map 
f(x) = rx(l — x). 

Following the analysis of the experimental data, the 
stochastic process of interest is represented by a coarse- 
grained variable defined as 

a n (t) = sgn[x n {t) - xtbx], t = 0, 2,4, . . . (2) 

where the quantity x n {t) is defined in Eq. Q and Xthr 
is a certain threshold value. The basic dynamics of the 
coupled oscillators being period two, the analysis is also 
done over even (or odd) time steps. 

This model shows a stretched exponential distribution 
of the renewal (or trap) times, which are defined in the 
next section III Bl over a wide range of the control pa- 
rameter r, which plays a role akin to temperature as 
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P to decrease with decreasing r (or T) , but also the quan- 
titative agreement in the low-r (or T) limit f3 « 1/3; the 
same limit as that obtained in studies of random walks 
on fractals [HEEl. 

Although the bifurcation diagram for a site in the cou- 
pled map lattice (Fig. O^d)) appears to be uniform with 
respect to r down to r rs 3.8275, the non-monotonicity 
of 0o (r) and To(r) most probably arises from the fact 
that for the isolated logistic map and the values of r < 4 
chaotic orbits are interspersed with period cascades seen 
in Fig. [He), i.e. the approach to full chaos as r goes to 4 
is not uniform. (For an explanation of the concept of bi- 
furcation diagrams see e.g. Ref. Interestingly, the 
lowest value of (3q is attained within the widest window 
of the period-three cascade in the bifurcation diagram of 
the uncoupled logistic map. About and below r = 3.82, 
in the coupled map chain, very narrow chaotic windows 
are also interspersed with periodic ones and the stretched 
exponential behavior of the trap-time distribution, with 
Po greater or about 1/3, is recovered only for some val- 
ues of r. This interesting behavior suggests that SER is 
uncovering a new aspect of the dynamics of coupled map 
lattices. A full analysis of this new behavior lies outside 
the framework of this paper, but we plan to return to 
studying it elsewhere. 

In other models and experiments, the time scale r 
monotonically grows with decreasing temperature. It is 
not the case here, although the values of To are greatest 
at the low-r (and low-/3o) end, following the general trend 
observed elsewhere. 



B. Renewal processes in coupled maps 



FIG. 1: /3q (a) and to (b) as functions of the logistic map 
parameter r. Dots correspond to the least squares-fit values 
and the thin solid line is drawn as a guide to the eye. The 
vertical dashed lines demarcate the discontinuities in /3o(r) 
and ro(r). They correspond very well to the positions of the 
periodic cascades in the bifurcation diagram of the uncoupled 
logistic map shown in (c). (d) shows the bifurcation diagram 
for a site in the coupled map lattice. Inset in (b) magnifies 
the right hand side of the main plot where the variation of 
ro with r is not distinguishable. The accuracy of /3o and to 
obtained by a least-squares fit to the long time part of the 
trap-time distribution is within 15%. 

can be seen in Fig. ^ This figure shows the stretch- 
ing exponent flo and the time scale To as a function of 
the control parameter r. The values of /3q and To were 
obtained by least-square fitting the long-time part of the 
trap-time distributions for a chain of N — 10, 000 cou- 
pled logistic maps f(x) — rx(l — x) with the expres- 
sion A exp[— (t/r )' 30 ]. Although /3 (r) andro(r) are non- 
monotonic, the overall behavior is similar to the results 
in both 3D [2j| and ID spin-glass models. This broad 
similarity not only captures the qualitative tendency for 



The stochastic process <J n (t) is a renewal process [23J 
provided that the time intervals (renewal times) between 
two subsequent zero crossings are statistically indepen- 
dent random variables. We tested the independence of 
these time intervals for the present model, Eq. J2J), by 
calculating the distributions of the time intervals follow- 
ing time intervals of a specified length. The distributions 
were identical regardless of the length of the preceding 
time intervals. 

In general, a renewal process is defined so that a(t) = 
+1 for t < t < ti, a(t) = -1 for h < t < t 2 , a(t) = +1 
for t 2 < t < £3, and so on. The renewal processes are the 
simplest possible stochastic processes readily susceptible 
of probability theoretical analysis (3^1 . The dynamics 
represented by such processes can be understood in terms 
of transition probabilities between the two possible states 
a = ±1 0. 

In [I] , it is shown that the existence of very long renewal 
time intervals t n+ \ — t n associated with the stretched ex- 
ponential distribution can be attributed to the presence 
of time-limited ordered domains (or traps) with domi- 
nant spatial period of 4 sites. Fig. OJa) shows a typical 
distribution of the renewal time intervals (or trap times) 
for this model. 
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FIG. 2: (a) Second derivative C"(t) of the autocorrelation 
function C{i) of the coarse-grained orbit a{t), Eq. J2) with 
r = 3.83, (solid line) compared to the asymptotic behavior 
p(t)/(t) given by Eq. J5J (dashed line). Dash-dotted lines 
shows the stretched exponenetial fit to the long-time part pf 
p(t) with T = 59±9 and (3 = 0.39±0.05). (b) The same as in 
(a), but for the renewal process computer generated with the 
renewal times distributed according to Eq. (JIJ with to = 5 
and /3o = 0.4. In this case, the lines corresponding to the 
dashed and dash-dotted lines in (a) coincide. 
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FIG. 3: (a) Autocorrelation function C(t) of the coarse- 
grained orbit a(t), Eq. @ wiht r = 3.83, (solid line) and 
the stretched exponenetial fit to the long-time part of C(t) 
with t c ~ 1250 and f3 c ~ 0.58 (dashed line). The inset shows 
the same curves, but with the logarithmic vertical axis, (b) 
The same as in (a), but for the renewal process computer gen- 
erated with the renewal times distributed according to Eq. J3J 
with to = 5 and (3q = 0.4. In this case the stretched exponen- 
tial fit to C(t) is with t c « 69 and /3 C » 0.54. 



These results suggest that the dynamical heterogene- 
ity manifested by the presence of time-limited ordered 
domains could be at the origin of the stretched ex- 
ponential relaxation in supercooled liquids and glasses 
and the renewal processes could provide a simple pic- 
ture for describing heterogeneous dynamics in these con- 
densed phases. In order to support these two affirma- 
tions, we first need to show that (1) there is a one to 
one correspondence between renewal processes and their 
time auto-correlation functions, because experimentally 
measurable relaxation responses are mathematically de- 
scribed by auto-correlation functions of certain dynami- 
cal variables and not directly by non-observable renewal 
processes, and that (2) the renewal processes are present 
in supercooled liquids and glasses and can be calculated 
effectively, exhibiting unambiguous stretched exponential 
distributions of renewal time intervals. 



The relation between the auto-correlation function 

C(t) =< a(t')a(t'+t) > v (3) 

of a renewal process a(t) and the distribution p(t) of 
time intervals t n+ i — t n between the zero crossings (or, 
in other words, renewals) of this process was recently 
studied for a range of distributions [2jJ. Following this 
approach, Godreche and Luck also showed [3_5J that C(t) 
corresponding to the stretched exponential distribution 

Pit) = exp[-(t/ro)^] (4) 

can also be represented by a stretched exponential at long 
time 

C(t) w C s . c .(t) = Aexp[-(t/r)% t -» oo (5) 

We can Laplace transform the distribution of inter- 
vals p(t), obtaining a function with a small singular part, 
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p(u) = 1 — (t)u + ••■ + p s i n g(u), u — > [36J, where 
(t) = r r(2//3o]/r(l//3 ), is the first moment of p{t). 
Also from Ref . |22( , we know that the Laplace transform 
of the auto-correlation function is given by 



Ceq(lt) 



1 



1 - 



2(1 -P{u)) 
{t)u(l + p(u)) 



(6) 



Using the expression above for p(u), one finds, for the 
singular part of the correlation function at equilibrium 



Ceq,sing (^) 



Psing(u) 

(t)u 2 



(7) 



Since the behavior at large time is governed by its singu- 
lar part, we get 



C"(t)^p(t)/(t). 



(8) 



This implies that for a stretched exponential p(t) in the 
asymptotic limit C eq (t) oc yo(t)[3^|. 

The range of relevance for this analytical result can be 
confirmed numerically by studying the statistical proper- 
ties of a computer-generated sequence of (pseudo)random 
numbers distribution according to the stretched ex- 
ponential probability density function (pdf) Eq. I@J). 
Using the fundamental transformation law of proba- 
bilities |37| . given a computer-generated sequence of 
(pseudo)random numbers {9} distributed uniformly in 
the interval < 9 < 1, we apply a transformation rule 
t = 7o[Q -1 (l//3o, 9)] 1 / 130 so that the resulting sequence 
{t} is distributed according to the stretched exponential 
pdf p(t). In order to compute the inverse regularized 
incomplete gamma function Q~ 1 (a,9) 1 we used the soft- 
ware from Refs. |38ll3^|. 

Given a sequence of time intervals distributed accord- 
ing to a prescribed pdf, it is straightforward to construct 
the corresponding renewal process and its autocorrelation 
function and second derivative. The second derivative is 
calculated numerically, after applying a smoothing func- 
tion before each derivative to improve the quality of the 
operation. 

The resulting C"(t) and C(t), shown in Figs. |21 and 
respectively, demonstrate clearly that the correspon- 
dence between p{t) and C"(t) holds for long times for 
C(t) calculated both for the coupled map lattice model, 
Eq. P with f{x) = rx(l - x), r = 3.83, a = 0.25 
and N — 1000, and for a renewal process where the re- 
newal time-interval distribution is stretched exponential 
by prescription. C"(t) obtained numerically agrees with 
p(t) I (t) , Eq. JSJ , without any free parameter within a 
multiplicative factor of the order of unity. 

Interestingly, for the range measured, C(t) also follows 
the stretched exponential form, albeit with somewhat 
different parameters r c and f3 c . This indicates that the 
true asymptotic regime for C(t) or p(t) may be reached 
at longer time with a gradual approach towards this 
long time behavior. This agrees with the theories (see 
e.g. [24|]) where the stretched exponential functions are 



accompanied by time-dependent multiplicative factors. 
For practical purposes, however, it is remarkable that 
the stretched exponential fit to the correlation function 
agrees with C(t) starting with a well-observed experimen- 
tally value between 3/4 and 1/2 x max{C(i)} and follows 
with very high accuracy for several orders of magnitude 
to the regime beyond the reach of the experiments. 



C. Molecular dynamics 

Having established the correspondence between the 
distribution of time intervals of renewal processes and 
their auto-correlation functions, we now show that a re- 
newal process can also be identified in a supercooled liq- 
uid and it appears to provide a simple and elegant de- 
scription for the stretched exponential relaxation. In 
order to do so, we have performed molecular dynam- 
ics simulations of a 16000-particle single-component sys- 
tem where the interactions between the particles are 
described by a pair potential labeled Z2 in Ref. [40| . 
This potential is similar to that introduced initially by 
Dzugutov but provides a longer-lived metastable su- 
percooled state prior to crystallization. Otherwise the 
properties of the Z2 liquid are similar to those of the 
Dzugutov liquid, in particular it exhibits dynamical het- 
erogeneities formed by short-lived clusters composed of 
connected icosahedra [42^. The simulations are per- 
formed at T = 0.68, below the melting temperature of 
T m = 0.70 ± 0.05 with a time step of 0.01 in reduced 
Lennard- Jones units We take measurement over 

the 16 000 atoms with the number density p = 0.85 for 
7 million time steps, stopping before the system starts 
crystallizing. This simulation time would correspond to 
a 150-ns simulation for argon. 

A renewal process representing the atomic dynamics 
in the liquid can be constructed so that each renewal 
time interval correspond to the time a particle spends in 
a dynamical trap. This is identical to the first passage 
time, in the parlance of Allegrini et al. 0], but differs 
formally from the waiting time of Doliwa and Heuer [l9| 
and Denny et al. [2(j which focuses on hops between en- 
ergy basins. We say that particle i is trapped as long as 
its displacement 



dj(t) = |r;(i) -r,(0)| 



(9) 



measured starting from an arbitrary initial position rj (0) 
remains less than a certain threshold value dthr- The 
moment t — i t hr of surpassing the threshold is counted 
as a renewal and the value of ri(0) is reset to rj(tthr) as 
well as the time origin for the next trap. 

The corresponding trap-time distribution (see Fig. 0J) 
is not particularly sensitive to the choice of dthr as long 
as it is of the order of the effective atomic diameter. For 
definitiveness, we choose d t hr = 2n/Q ~ 0.88 ^3|> where 
Qo is the position of the main peak in the static structure 
factor S(Q) and neglect the temperature dependence of 
this length scale. A threshold of « 0.88 is about 4 times 
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larger than the critical value estimated by Allegrini et al. 
in a binary Lennard- Jones system. While the focus of 
Ref. ^| is on the inertial regime, which takes place as 
small distance, on the order of 0.1 atomic radius, we are 
interested in the atomic displacement leading to a change 
in the configuration. Figure 0{a) also shows in inset that 
the trap-time distribution is exponential at the temper- 
ature well above the melting point (T m — 0.70 ± 0.05) 
and its long-time tail well agrees with a stretched expo- 
nential fit with Pq w 0.71 in the supercooled regime. In 
order to be considered as a fully Markovian process, we 
also assess the correlation between subsequent trapping 
time intervals, as was done for the coupled map lattice. 
Comparing trap distributions considering only the subset 
of events following traps of various length, again, we find 
no correlation between renewal time intervals. The value 
of Pq agrees well with j3 = 0.70 ± 0.05 for the self part 
of the intermediate scattering function F s (Qo,t) calcu- 
lated at the same temperature and for the same length 
scale. (For the definition of F s (Q,t) see, e.g., Ref. _44j.) 
Fs(Qo, t) for the Z2 liquid behaves similarly to the same 
quantity calculated for the supercooled Dzugutov liquid 
l45l and is not shown here. 



These results appear at odds with the analysis of Al- 
legrini et al. who report a power-law distribution for 
the long-time tail of the renewal time i although ad- 
mitting that the asymptotic cutoff in the distributions 
is difficult to resolve numerically. Fig. Bib) shows our 
simulation results at T = 0.68 plotted on a log-log scale 
with a tentative fit of the long-time behavior. It is clear 
that the stretched exponential provides a much better 
fit to the simulation data. This difference between the 
two simulations can be due to both a physical range of 
parameters and the precision of the simulation. First, Al- 
legrini and collaborators focus on the inertial length-scale 
while the results presented here require atomic diffusion, 
not only thermal vibrations; second, focusing on a single 
temperature, our simulation uses a larger cell and aver- 
ages over a time scale roughly 6 times longer. 
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FIG. 4: (a) Dots: distributions of trap-time intervals for 
low (T=0.68, each data point corresponds to an average 
over 500 time steps) and high temperatures (inset, T=2). 
Dashed lines: corresponding stretched exponential fits with 
r = 60 ± 5, fa = 0.71 ± 0.05 and (inset) r = 1.22 ± 0.05, 
Po = 1 ± 0.002. (b) Dots: the same as in (a) for T=0.68 pre- 
sented on a log-log scale over a mesh with logarithmic time 
intervals. Dashed line: the same as in (a) for T=0.68. Dash- 
dotted line: power law tx 1/t . 



III. CONCLUSION AND DISCUSSION 



Denny et al. |20j also note that the distribution in 
Ref. hj may not be a power law. Although the wait- 
ing time distribution in Ref. |20j | differs formally from 
that is studied here, at long time they can be compara- 
ble since the escape of a particle from a sphere of radius 
of the order of the atomic diameter is likely to correspond 
to the transition of the phase point in the configuration 
space between energy metabasins. We note however, that 
the log-normal form of the waiting time distribution pro- 
posed by Denny et al. results in a correlation function 
that can only be approximated with SER, whereas the 
stretched exponential first passage-time distribution re- 
sults in a correlation function that is also stretched ex- 
ponential at long time. 



As was shown in [I] , traps in coupled map lattices corre- 
spond to ordered short-lived domains. In direct analogy, 
it is natural to expect that the trap condition di(t) < dthr, 
for di(t) from Eq. is satisfied as long as particle i is 
literally trapped within a relatively stable environment 
- a "cage". Thus formed renewal process and the as- 
sociated trap-time distribution and the autocorrelation 
function provide a simpler description of the structural 
relaxation in a condensed matter system than, e.g., the 
cage correlation function recently used to demon- 
strate the stretched exponential relaxation in a Lennard- 
Jones massively defective crystal. 

In particular, traps in the coupled map lattice show 
only very short spatial extension that does not seem to 
diverge as the overall dynamics of the system slows down. 
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As (3 goes from 0.6 > (3 > 0.5 to 0.4 > /3 > 0.3, the 
longer traps go from about 120 to more than 30,000 time 
steps but the characteristic time-averaged width of these 
traps only quadruple, from 2 to 8 sites [I]. These ob- 
servations are consistent with the behavior of dynamical 
heterogeneities observed in a number of model systems 

We note that the only model that rigorously derives 
the stretched exponential long time asymptotic behav- 
ior is that of a random walker in a system with static 
traps (see t^g. Ref. 0). Ref. cites both upper 
and lower 48] bounds for the asymptotic behavior of the 
survival probability. Both bounds have the same /3, but 
different r. It is not obvious, however, that the above 
theory directly applies to the results of this works. Nev- 
ertheless, our results are obtained in terms of first pas- 
sage statistical distributions of the type encountered in 
the static-trap model and the use of similar mathemati- 
cal techniques may shed more light on the nature of the 
dynamics in the models presented here. 

In conclusion, we focused on a relation between dynam- 
ical heterogeneities and stretched exponential relaxation 
in two models: a coupled chaotic oscillator system and 
a supercooled liquid. Without insisting that the former 
model represents the actual behavior of real liquids, we 
have shown that the asymptotic long time behavior of the 



autocorrelation function of a renewal process is stretched 
exponential provided that the distribution of the renewal 
time intervals p(t) is also stretched exponential at long 
time, i.e. p(t) ss exp[— (t/ro)' 30 ]. Using this relation, 
we have demonstrated that, at least in the two model 
systems considered here, the relaxation dynamics can be 
described in terms of well-defined dynamical traps, pro- 
viding useful insight regarding the absence of a diverging 
static length scale at the glass transition. 
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